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We present a new approach to the calculation of thermodynamic functions 
for crossing-invariant models solvable by Bethe Ansatz. In the case of the XXZ 
Heisemberg chain we derive, for arbitrary values of the anysotropy, a single non- 
^| linear integral equation from which the free energy can be exactly calculated. The 



high-temperature expansion follows in a sistematic and relatively simple way. For 
low temperatures we obtain the correct central charge and predict the analytic 
structure of the full expansion around T = 0. Furthermore, we derive a single 
non-linear integral equation describing the finite-size ground-state energy of the 
Sine-Gordon quantum field theory. 
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The computation of thermodynamical functions for integrable models started 
with the seminal works of C.N. Yang and CP. Yang [1], of M. Takahashi [2] and 
M. Gaudin [3]. In Refs. [2,3] the free energy of the Heisenberg spin chain is written 
in terms of the solution of an infinite set of coupled nonlinear integral equations, 
derived on the basis of the so-called "string hypothesis" . In this note we propose a 
simpler way to solve the thermodynamics by means of a single, rigorously derived, 
nonlinear integral equation. We restrict here our attention to the XXZ spin chain 
but the method can be applied to a wide class of models solvable by Bethe Ansatz. 

The XXZ hamiltonian for a periodic chain with 2L sites (we assume periodic 
boundary conditions but generalization to other b.c. is possible) takes the form 



2L 

HXXZ = -JY1 [«+! + ^l+l - COS 7 + 1)] (1) 

n=l 

The anisotropy parameter 7 is assumed here to lay in the real interval (0, ir) 
characteristic of the gapless regime. Hxxz is related to the diagonal-to-diagonal 
transfer matrix Ti(6) of the six- vertex model by 

T L (0) 6 =° 1 - tt^—Hxxz + 0{6 2 ) (2) 
2 J sin 7 



where 



Tl{9) = -R12-R34 • • • R2L-1 2L-R23-R45 • • • R2L 1 (3) 



and 



a + c a c j_^(-x„x ,^y„y \ 

Rum =— ^ 1 On°m + T^n^n+l + ^n^n+l) 



sm(7 — 6) j sin 

a= , b= , c=l 

sin 7 sin 7 



(4) 



9 is known as spectral parameter while a, b, and c are the conventional six-vertex 
Boltzmann weights. 
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The free energy per site is defined as usual 



fW) = ~- L lim -L log [Tr (e-^«)] (5) 



/3 L^cx) 2L 

and from eq. (2) we read 



e -f3H xxz = Um 



N 

T L (2(3/N) 



where (3 = flJsmj. Hence we can write 

f{(3) = _1 Um J_ lim logZ LiV (2^/iV) (6) 

p L^oo ZL N^oo 

where Zln(@) = T^[Tl(9)] N is the six-vertex partition function on a periodic 
diagonal lattice with Lx N sites. The two limits in eq. (6) cannot be interchanged 
since the degeneracy of T^(0) = 1, that is 2 2L , is strongly L— dependent. However, 
the crossing invariance of the six-vertex R— matrix implies that under a rotation by 
7r/2 of the entire lattice plus the substitution # ^ 7 — the numerical value of the 
partition function does not change (see e.g. ref.4). Therefore Zln(@) — Znl(j — 9) 
and 

f({3) = -- lim -L Jim logZ jVi ( 7 - 2/3/7V) (7) 

p L^oo ZL iV— >oo 

where Z^l(9) = Tr [T/v(#)] L . Now, as the well-known Bethe Ansatz solution tells 
us, Tjv(7) has a non-degenerate largest eigenvalue for any finite N. Then the two 
limits in eq. (7) commute and one finds [5] 

-2(3f((3)= lim logA£r( 7 -2p7A0 (8) 

where K^ ax (9) denotes the largest eigenvalue of Tn(9). Using the results of ref.[6] 
for the eigenvalues of the diagonal-to- diagonal tranfer matrix, the free energy 
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becomes 



><0 = ^J?5o 



E N (P)-2Nlog 



sin(2/?/iV) 
sin 7 



where 



(9) 



N 

e(A, (5) = <t>{\j + 2 7 /2, 7/2 - P/N) - <f>(\j - i 7 /2, 7 /2 - /3/7V) 

/ / ■> \ sinh(ix + A) 

0(A, x) = z log — — r- (0(0, x =0 

smh(zx — A) 



(10) 



The real numbers Ai, . . . , Ajv are the roots of the Bethe Ansatz equations 



N 



N 



<j>(\jJ/N) + d>(\j, 7 - P/N)\ = <M A i - X k, 7) + tt(2j — N — 1) (11) 

fc=i 



associated to the largest eigenvalue Ajy a:E ( 7 — 2(3 /N) (we assumed iV even). We 
define as usual the counting function [7] 



z N (X) 



2tt 



1 N 

0(A, P/N) + 0(A, 7 - P/N) ~j^J2^ X - X ^ l) 



fe=i 



(12) 



which by construction enjoys the properties zn(X) = —zn(—X) = zn{X) and 



z N (\j)=j-(N + l)/2 



(13) 



Our aim is to derive an integral equation for zn(X). 

Let us start by assuming that the root distribution becomes the continuus 
density cr c (A) (normalized to 1) in the limit N — > oo. Then it could be identified 
with the derivative of Zjv(A) and at the same time used to replace the discrete sum 
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in eq. (12) with an integration. For large but finite N we whould therefore obtain 
the following approximate linear integral equation for a c (X) 



+00 



2na c (X) = 0'(A, P/N) + <f>'(\, 7 - 0/N) - J ' dfi <f>'(\ - fi, 7 )a c (/i) (14) 



Fourier transforming leads to the explicit solution 



<7 C (A) = z' c (X) , z c (X) = — arctan 

7T 



sinh(7rA/7) 



(15) 



.sin(7r (3/>yN)_ 

We see that <x c (A) is strongly peaked at A = for large N, reflecting the singularity 
which develops at the origin in the iV — > 00 limit of the source term in eqs. (11), 
(14). The density picture correctly describes the BA roots wherever <r c (A) is of 
order 1. From eq. (15) we then find as validity interval 



|A| <0(v#V) 



(16) 



which shrinks to zero when N — > 00. Roots |Aj| > 0(y/~/3/N) have a spacing of 
order larger than 0(1 /N) and cannot be descried by densities. In particular, the 
roots with largest magnitudes have finite iV — > 00 limits spread by 0(1) intervals 
(we checked this fact numerically too). Therefore, contrary to the usual situation 
[7], we must go beyond the density description to obtain a bulk quantity like the 
free energy per site. 

Through simple manipulations we find from eqs. (12) and (14) [8] 



+00 



z' N (X) = z' c (X) - J dfipiX-fi) 



1 N 

jf^2S(n-Xj)-^ N (fi) 



3=1 



(17) 



where 



+00 



P(A) 



dk 
2^ 



p(k)e ikX , p(k) 



sinh(7r/2 — j)k 



(18) 



2sinh(vr - rfk/2 cosh(7fc/2) 

—00 

The r.h.s. of eq. (17) can be expressed as a countour integral encircling the real 
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axis, that is 

+00 



z'n( X ) = 4(A) + J dfip(\ - fi)z' N (fi) 



+ 



\ _|_ e 2niNzN(ft—iO) J _|_ g— 2niNzN(n+iO) 

(19) 

where we used eq. (13). Integrating eq. (19) now yields 

+00 

log— — = dfip(X-fi) log- — (20) 

a c (A) y l + ajv(A* + iO) 

—00 

which is the sought nonlinear integral equation written in terms of ax = exp(2iriNz jv)| 
(similarly a c = exp(2niNz c )). Notice that ajv(A) = l/ajv(A). We evaluate the sum 
in eq. (10) by a similar procedure, with the result 

+00 

r . . ^ 1 „ . _ If,, sinh7rA/7 costtBHN , l + a/v(A + iO) 

L N (P)= E N (P) - E C (P) = - dX ^ ^4 log v ; 

7 7 cosh 2ttA/ 7 - cos 27r/5/ 7 AT 6 1 + a N (X + iO) 

—00 

(21) 

where E c {j3) is that part of the energy that follows from the root density (15), that 
is 

+00 +00 
^/^x f x /" ^ sinh(7r — 7) A; sinh(7 — 2d/N)k ,_ x 

— OO —OO 

If eqs. (20) and (21) are analytically continued to the axis ImA = ±7/2 they 
assume the same structure of those derived in ref. [10] by different methods. 

Inserting eqs. (21) and (22) in eq. (9) yields 

f(P) = E xxz ( 1 )+p- 1 m (23) 
where Exxz{l) is the ground-state energy of (1), namely 



Exxzil) = 2J 



+00 



/ sinh(7r — 7) A; 
dk 



sinhvrA; cosh 7 A; 

OO 



(24) 
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while L(/3) is the N — > oo limit of (21), that is 



+00 



y 7sinhvr(A + iO)/7 v ; 

— oo 

Similarly, the function a(A) in this expression is the iV — > oo limit of ajv(A). It is 
therefore the solution of the N — > oo limit of eq. (20), that is 



+oo 



-iloga(A) = ~~~~~ . + 2 / dup(A - u) Im logfl + a(u + iO)l (26) 

7smh7rA/7 J 

— oo 

Thus the calculation of the free energy (5) has been reduced to the problem of 
solving the single nonlinear integral equation (26) and then evaluating the integral 
(25). 

It is possible to rewrite eq. (26) in the alternative form 

+oo 

z(X) = (3q(X) + ^2 Im J dfj,<j>'(\ -ji,7)logcos7rz(ji + i0) (27) 

— oo 

where 

1 ( sinh2A , A 

q(X) = - — — - coth A 28 

Hy ' it Vcosh2A - cos2 7 J v 1 

and z(X) = 2^1oga(A) is related to the original counting function 

z(X) = lim N \zn(X) - isign A] (A real) (29) 

By the residue theorem we then obtain 

z(X) = 0q(\) - i- m - 0,7) - 0(A, 7)] (30) 

where the real numbers £j, defined by z(£j) = j — 1/2, j e Zi, can be identified 
with the N oo limit of the original BA roots Ai, A2, . . . , Ajv- Eq. (30) shows that 
z(X) has periodicity in and has, as unique singularity on the real axis, a simple 
pole at the origin with residue 
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Let us now study f(/3) for high temperatures. When f3 is small it is convenient 
to use the alternative form (27) plus the uniform expansion 

oo 

z{\) = Y,pb h {\) (31) 
fc=l 

Since the residue of z(X) is linear in (3 and z{6) is an odd function, we have 

6 1 (A) A ~°-^ + 0(A), M0) = (k>2) (32) 
Inserting eq. (31) in eq. (27) we find, up to fourth order in (3 



6i(A) = g(A) & 2 (A) = -V'(A, 7 ) 



(33) 



MA) = MA)^Q-^)^(A,7) + ^(A,7) 

Then, from eqs. (23) and (25), 

= _/r 1 ln2 + J cos 7-/? J 2 (l + ±cos 2 7) + /3 2 J 3 cos 7 + 0(/3 3 ) (34) 

which indeed agrees with the high T expansion [2]. We want to remark that eq. 
(27) generates the functions 6fc(A) recursively, with easy integrations which involve 
only delta functions and derivatives thereof. It is indeed a very efhcent way to 
recover the high-temperature expansion from the Bethe Ansatz solution [11]. 

We shall consider now the low temperature regime. When (3 ^> 1 eq. (26) 
indicates that z(X) ~ (3, so that log[l + a(A)] ~ 2mz(X), at least as long as 
|A| < (7/77) log p. At dominant (3 3> 1 order eq. (26) linearizes in the same way as it 
does for small (3. Inserting then z(X) ~ (3q(X) in eq. (25), yields lim^oo L((3) = 0. 
Therefore eq. (23) tells us that 

f{(3) P =°° Exxz^ + Oir 2 ) (35) 

The contributions of order (3~ 2 and smaller come from values of A larger than 
(7/71") log (3, where the previous assumption log a ~ 0((3) does not hold anymore. 
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It is then convenient to introduce the new function 



7 47T/3 

A(x)=a(X), x = \--\og — - 

it 7 



Then eqs. (25) and (26) reduce, in the f3 — > oo limit, to 



+00 



L((3) = -L [ dx e~" /7 Im log[l + A(x + iO)] 
ir/3 J 



(36) 



(37) 



and 



+00 



-ilogA(x) = -e wx ^ + 2 J dyp(x-y)lm\og[l+A(y + iO)} (38) 

— oo 

Fortunately, it is not necessary to solve this integral equation to calculate the 
integral (37). In fact, we can use the following lemma. 

Lemma. Assume that F(x) satisfies the nonlinear integral equation 



X 2 

—i log F{x) = ip(x) + J dyp(x — y) Im log[l + F(y + ii 



(39) 



where (p(x) is real for real x and x\, xi are real numbers. Then the following 
relation holds 



Im J dx<p'(x) log[l + F(x + iO)] = ilm [<p(x 2 ) log(l + F 2 ) - <p{xi) log(l + Fi)] 

Xl 

+ iRe [m) - £(F 2 )} 
where F\^ = F(xi j2 ) and £ is a dilogarithm function 



(40) 



l{t) = J du 



loerfl + u) logu 



u 



1 + u 



(41) 
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To prove the lemma we consider the relation 



x 2 

e{F 2 )-£{F 1 ) = J dx^og[l + F{x)]^\ogF{x) - \ogF{x)^\og[l + F{x)] 



XI 



and then use eq. (39) and its derivative to substitute logF(x) and d log F(x)/dx. 
After taking real part and a little algebra this yields eq. (40) (related identities 
were used in ref. [10]). 

The integral in eq. (37) may now be exactly calculated by invoking the lemma 
with F(x) = A(x), ip = — exp(— irx/^f) and x\ = — oo, xi = +oo. We have 
A(xi) = 0, A(x 2 ) = 1 and 

+oo 

2 Im I dx e~™^ log[l + A(x + i0)} = --£(1) = - — 

J 7T 6 

— oo 

Then the free energy for low temperature reads 

f{P) = E X xz{l) - Trr—r 2 + o(/T 2 ) (42) 
6 J sin 7 

in perfect agreement with refs. [2,9]. As we have just shown, both the high and the 
low temperature leading behaviors of the free energy can be derived without much 
effort from our non-linear integral equation (26). The higher order corrections for 
high T can be obtained in a very sistematic way [11]. The situation for the o((3~ 2 ) 
terms in the low T expansion (42) is more involved. Qualitatively, however, it is 
rather easy to establish, from the asymptotic behaviour of the kernel p(X — ji) in 
eq. (26) and from the in periodicity implied by eq. (30), that these higher order 
terms must be integer powers of T 71 " 7 and T 7 ^ 71 " -7 ). 

Let us now consider the problem of calculating the (properly subtracted) 
ground state energy E(L) of 2D integrable massive models of Quantum Field 
Theory in a finite 1-volume L. This subject recently received much attention 
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in connection with Perturbed Conformal Field Theory. In this contest E(L) is 
known as ground-state scaling function. Our starting point is the lattice regular- 
ization of such models provided by the light-cone approach to integrable vertex 
models [6]. Specializing to the six-vertex model, the ground state energy, on a ring 
of length L formed by N sites, takes a form similar to eq. (21): 



+00 

E N {L)-E C {L) = -^Im j d\ 

—00 



sech ^(A - 0) - sech 2(A + 6) \og[l+a N (\+iQi)} 



(43) 



where 



N 2 

Ec(L) = X 



+00 
-2tt+ J d\ 

—00 



4>(\ + 20,7/2) 



7 cosh 7rA/7 

and ajv(A + iO) obeys an equation like (20) with the new source term 



(44) 



loga c (A) = 2arctan 



sinh 7rA/7 
cosli7r6/7 



(45) 



The L— dependence in these expressions is hidden in the light-cone parameter 6, 
which tends to infinity, in the continuum limit iV — > 00, as [6] 



e = - log — 

7r mL 



(46) 



with m, the physical mass scale, held fixed. Define now 



e(A) = — lim loga7v(7A/7r) 

N^oo 



(47) 



then from eqs. (20), (42-45) we find 



+00 



E(L) = lim [E N (L) - E C {L)\ = -m I —sinh Aim log 

N^oo J IX 



1 + e -e(A+iO) 



(48) 
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and 

+00 

e(A) = mL sinh A + 2 J dfiG(X-fi) Im log 1 + e^+^j (49) 

—00 

where G(X) = ('j/n^^X/n). Thanks to the correspondence between the light- 
cone six-vertex model and the sine-Gordon (or Massive Thirring) model [6], eqs. 
(48) and (49) give the ground state scaling function E(L) of the sine-Gordon field 
theory on a ring of length L. This is a rigorous and simpler alternative to the 
standard Thermodynamic Bethe Ansatz [13]. Notice that all UV divergent terms 
are contained in E C (L) (eft. eq. (44)). E C (L) also contains the finite, scaling part 
m 2 Lcot 2 7r 2 /27 (also found in ref. [12] by completely different means), which gives 
the ground-state energy density in the L — > 00 limit. 

Let us now study eqs. (48) and (49) for mL 3> 1 and mL <C I. For small 
values of mL the calculation closely parallels that for low T in the thermodynamics 
of the XXZ chain. Shifting A — > A — log mL and applying the lemma leads to 

£(L) m 4r°--^ + 0(l) (50) 

showing the expected value c = 1 for the UV central charge. The large mL regime 
easily follows from eq. (49) by iteration. One finds the behavior typical of a massive 
quantum field theory 

E(L) L ~°° -— Kl (mL) + 0(e- 2mL ) (51) 

7T 

with Ki(z) the modified Bessel function of order 1. 

We would like to remark that, contrary to the traditional thermal Bethe Ansatz 
[2,3], our approach does not relay on the string hypothesis on the structure of 
the solutions of the BA equations characteristic of the Heisemberg chain. This 
makes our approach definetly simpler. Most notably, the whole construction of 
the thermodynamics no longer depends on whether 7/7? is a rational or not, unlike 
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Takahashi approach [2]. This applys equally well to the problem of the ground 
state scaling function of the sine-Gordon field theory, since the standard TBA 
approach requires the string hypothesis. 

In this note we restricted ourselves to the XXZ Heisenberg chain and to the 
sine-Gordon model. Generalizations of our approach to higher spin chains as well 
as to magnetic chains and quantum field theories associated to Lie algebras other 
than A\ (nested BA solutions) are relatively straightforward provided crossing 
symmetry holds [11]. 
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